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Abstract 


Music can be viewed as a structure formed by notes. Different structures in music have potential to yield 
enormously rich diversity of different melodies. Music is a typical example where a structure defines 
a property. Similar concepts could be also exploited in education, in presented case in nano-sciences, 
which are typical representatives of soft materials the structure of which can be strongly manipulated 
with local geometry and presence of appropriate nanoparticles. The objects of study, named also LC 
shells, will be exploited as basic unit elements for future soft colloidal crystals. A different arrangement 
of colloids within the crystal would result in different physical properties in a similar way as different 
packing of atoms results in different crystals made of real atoms. 

In presented research will be demonstrated, how relevant basic mechanisms in thin films of nematic liquid 
crystals could be explained in a classroom and used as a case study, also for explanation of many other 
Physical properties. This research topic is still in its infancy. At this stage only various defect structures 
in relatively simple geometries (spherical and elliptical) will be analyzed. There is a need to find simple 
ways to control sensitively the valence of LC shells and in particular to develop strategies to assemble 
them in crystal structures of desired symmetry. This would allow tailoring specific optical dispersion 
relations or other physical property of interest and make new ways to teach different physical properties 
on the »music« based approach. 

Key words: teaching/learning strategies, »music« based approach, topological defects, localized states, 
nanoscience, soft materials, liquid crystalline shells. 


Introduction 


History manifests strong interplay between art and science. Roughly 2.500 years ago 
Pitagora was the first to observe connection between mathematics and music. Later Johann Se- 
bastian Bach discovered universality in music, which is nowadays broadly applied in science. 
Furthermore, in his several compositions he combined fine balance between “order” and “dis- 
order” which enabled him to obtain rich diversity in music expressions. Note that the key in- 
gredient of modern complexity theory, which explains richness of natural phenomena because 
the Nature is at the “edge of chaos”, separating regions exhibiting chaotic behaviour and order. 
Cubism in art, in particular Georges Braque, inspired pioneering research in Quantum Mechan- 
ics to introduce the duality concept of Nature, where either wave or particle- type behaviour is 
observed. At the beginning of 20" century Henri Poincaré introduced phase space in order to 
solve differential equations graphically. His mathematical three-body behaviour analysis trig- 
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gered discovery of Chaos, which later nucleated science of complex systems. Paintings of Ja- 
kson Pollock represent first artistic manifestations of fractals before they were discovered in 
science. After Chaos was accepted as an important branch in science, different fractal textures 
become ubiquitous and fashionable ingredient in paintings. Later Gyorgy Ligeti introduced 
fractals, representing the key signature of a chaotic behaviour, into his music compositions... 
These few key examples evidently show mutual connection between science and art, where 
experiences gained in one branch could be beneficial to the other branch. 

Physical science, particular nanoscience could also benefit from music. Namely, various 
melodies are in fact different structures composed of finite number of discrete units - notes. By 
changing their relative position, duration, dynamics (1.e. rhythm) different melodies and music 
styles emerge. This paradigm, where different structures define qualitatively different global 
behaviour, can be applied in nanoscience. Currently, main focus in constructing nano-based 
new materials is based on developing new nano-units. However, music suggests that different 
»»properties« are not always directly defined by the basic building blocks themselves but by the 
way they are organized into hierarchies. Music is in fact a hierarchical system, where specific 
patterns are nested within larger patterns. Similar mechanisms might be applied in designing 
new nano-composites with exotic properties. 

As an illustration of »music« based approach we present different structures in thin liquid 
crystalline (LC) films (De Gennes, Prost, 1993) and LC shells (Vitelli, Nelson, 2006). The latter 
represent thin films of nematic LCs covering a colloidal particle. As “units” we consider topo- 
logical defects (TDs) in LC orientational ordering (Mermin, 1976). Our research is theoretical 
and we use two dimensional Landau-type approach in terms of complex order parameter (De 
Gennes, Prost, 1993, Kralj, et al., 2011). Using it we express the free energy density in terms of 
invariants formed by the Taylor expansion in Q and its first spatial derivatives. Different metast- 
able or stable structures of orienational ordering defined by Q are obtained numerically by 
solving corresponding Euler-Lagrange equations (Kralj, et al., 2011). Different patterns of TDs 
by exploiting their coupling with the curvature of LC films or shells and TDs are imposed and 
shown how such systems could be in future exploited to form the so called scaled crystals (Nel- 
son, 2002), where microscopic LC-shells are analogues of atoms in conventional crystals. 


Localised States in Continuum Theories 


In nature localised “objects” with well defined fixed shapes (1.e., constant in time) re- 
minding us of notes in music have been observed and referred as particles. At a fundamental 
level, in physics ordering of a system is commonly described in terms of continuum fields. 
Scientists have for years searched for nonlinear field-type theories yielding stable localised 
states, reminding us of particles. Skyrme was the first to describe nuclear particles as local- 
ised states (Skyrme, 1961), so called “Skyrmions”’. In the following years they discovered 
Skyrmions in field-type theories describing variety of completely different physical systems, 
where characteristic lengths varied from microscopic (Al’ Khawaja, Stoof, 2001) to cosmologi- 
cal scales (Durrer, et al., 2002). In these pioneering studies Skyrmions have been found only 
under non-equilibrium conditions or have been stabilised by external fields. However, latter 
stable configurations were discovered, for example in magnetic systems (ROBler1, et al., 2006), 
where Skyrmions are stabilised by chirality. Chirality refers to symmetry breaking between left 
and right handedness. Therefore, a chiral system 1s characterised by a preferred handedness and 
consequently the system lacks the inversion symmetry. 

In a wider context Skyrmions are members of the family of topological defects (TDs) 
(Mermin, 1976). Namely, their local stability is due to topological reasons. Consequently, their 
main property is defined by a topological quantity, so called topological charge, which 1s in- 
dependent of microscopic details of a system. TDs appear at least temporary in systems which 
are reached via a symmetry breaking phase transition. In fact, last majority of structures in 
condensed matter and even structure of the universe were reached via such transitions. For this 


ISSN 1822-7864 


Samo KRALJ, Boris ABERSEK, Irena KRALJ. Interplay between Art and Science in Education: “Music” Based Approach in Nano- 
sciences 


reason TDs are ubiquitous in nature, and science of topological defects is universal. For exam- 
ple, the first theory of coarsening dynamics of topological defects was developed in cosmology 
(Zurek, 1976). However, in most commercial applications based on TDs liquid crystals are used 
as working materials (De Gennes, Prost, 1993). 


Topological Defects in Liquid Crystals 


Skyrmion-type stable configurations are observed also in liquid crystals (LCs) (De 
Gennes, Prost, 1993). Most famous are Blue phases (BPs) (Cordoyiannis, et al., 2013) and 
Twist Grain Boundary (TGB) (Wright, Mermin, 1989) phases which are dominated by stable 
lattices of such topological defects. These are in BPs and TGBs commonly referred to as line 
disclinations and screw dislocations, respectively. Note that the discovery of the latter structure 
was guided by mathematical similarity in describing translational ordering in LCs and super- 
conductivity in superconductors (Chaikin, Lubensky, 1995). Namely, in both systems the pri- 
mary order parameter field is a complex quantityy =he" , where the modulus h (degree of 
ordering) and the phase factor f are real numbers. In this analogy 1) TGB and the 11) screw 
dislocation in LCs are analogues of the 1) Abrikosov phase and the 11) magnetic vortex in super- 
conductors. 

Liquid crystals (De Gennes, Prost, 1993) are typical soft matter representatives, dis- 
playing ordering between conventional liquids and solid crystals. The defining property of soft 
materials is easiness of their response to even modest external perturbations. Therefore, LCs 
can be relatively easily manipulated by confining surfaces and/or external electric or magnetic 
fields. LCs are optically transparent and anisotropic and consequently their structural character- 
istics could be relatively easily monitored using, for example, optical polarization microscopy. 
Furthermore, richness of different LC phases and structures enables existence of any topo- 
logical defects encountered in Nature (Mermin, 1976). Consequently, LCs are most convenient 
testing ground systems to study fundamental properties of topological defects. 

The nematic phase is the simplest LC phase (De Gennes, Prost, 1993). For presentation 
purposes we henceforth restrict to thermotropic LCs 1n which the nematic phase is reached from 
the isotropic (ordinary liquid) phase by lowering temperature. We also set that on average LC 
molecules behave as rod-like objects with head-to-tail symmetry (1.e., up and down positions 
are equivalent). In bulk samples (relatively large sample where we neglect boundaries) in the 
nematic phase LC molecules are on average aligned along a single (the so called symmetry 
breaking) direction. The nematic local average orientational ordering at a site 7 1s commonly 
described by the unit vector (7), called the nematic director field (De Gennes, Prost, 1993). 
The degree of ordering about this direction is quantified by the scalar nematic order parameter 
field | (7) (De Gennes, Prost, 1993). For rigidly aligned LC molecules it has maximal value. 
On the other hand | =0 reflects absence of orientational ordering, where all orientations are 
equivalent. 

In a nematic LC phase topological defects can be relatively easily visualized. Some 
simple examples are sketched in Figure 1. They refer to points or lines, where the director field 
is not uniquely defined. In most cases a large number of TDs appear temporary in relatively fast 
LC structural or phase changes enforcing local frustrations, which could not be relieved by 
continuous spatial variations inn(7). Note that local free energy costs of introducing a topo- 
logical defect are in general relatively large. Therefore, in most cases systems tend to avoid 
presence of TDs. However, due to topological reasons one could not get rid of an isolated TD if 
the director field on any surface enclosing the defect is fixed. One can decrease number of TDs 
only if a pair of defects of opposite topological charge annihilate each other (Kralj, et al., 
2011). 

Principal property of a topological defect is its topological charge m (Mermin, 1976), 
which is a topological invariant. Therefore, it does not change if we make continuous changes 
within a volume containing the defect if conditions at the volume’s boundaries are unaltered. In 


ISSN 1822-7864 


PROBLEMS 
OF EDUCATION 


IN THE 21 CENTURY 
Volume 61, 2014 





PROBLEMS 
OF EDUCATION 


IN THE 21* CENTURY 
Volume 61, 2014 





Samo KRALJ, Boris ABERSEK, Irena KRALJ. Interplay between Art and Science in Education: “Music” Based Approach in Nano- 
sciences 


our study we will restrict to two dimensional systems (2D). In 2D the topological charge is 
equivalent to the winding number (also referred to as the Frank index (De Gennes, Prost, 1993)), 
which can be simply visualized. The winding number is determined by the orientational field 
surrounding TD. It counts total angular change in orientation of 7 on encircling TD clock-wise 
using an arbitrary path divided by 20 . If TD exists then m equals to half integer or integer. 
TDs exhibiting m > QO and m <0 are commonly referred to as defects and antidefects, respec- 
tively. In Figure 1 are depicted schematic plots of simple TD structures bearing charges m = 1/2 
, m=-1/2, m=1 and m=-1. Note that TDs in most cases behave in a similar way as 
electrically charged particles or fundamental particles in particle physics bearing some quan- 
tized charge in general (ROBler1, et al., 2006). Namely, the total charge of a closed system is a 
conserved quantity. TDs of opposite (same) charges attract (repel) each other and can annihilate 
into a defectless state. 
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Figure 1: Topological defects in orientational ordering are presented by topo- 
logical charge which is in two dimensional systems equivalent to the 
winding number m. a) m=1/2, b) m=-1/2, c) m=1, d) m=-1. 


Soft Nano-Composites: “Music” Based Approach 


Recent years face increased interest into soft nano-composites (Nelson, 2002, 
Cordoyiannis, et al., 2013) and also on impact of nanoparticles (NPs) on biological membranes 
(Vitelli, Turner, 2004). Soft nano-composites typically consist of a soft matrix containing NPs. 
In such a way unique properties of soft materials and nanoparticles are combined. It is of interest 
that the resulting material either exhibits anomalously enhanced quantitative properties (with 
respect to “pure samples” in absence of nanoparticles) or qualitatively new properties. 

In particular it has been shown, that TDs within soft matrices could strongly interact with 
nanoparticles (Cordoyiannis, et al., 2013). Several studies reveal that appropriately surface 
decorated NPs could be efficiently trapped within cores of topological defects in various LC 
phases. Furthermore, appropriate NPs could substantially increase (by orders of magnitude) the 
stability range of relevant defect lattices due to the universal Defect Core Replacement (DCR) 
mechanism (Cordoyiannis, et al., 2013). This describes reduction of condensation free energy 
penalty related to existence of defect’s core by partially replacing its costly core by volume of 
NPs. This mechanism is efficient if NPs do not significantly disrupt a local LC environment. In 
such a way one could control positioning of NPs via manipulation of positions of TDs. The 
latter could be relatively easily manipulated in different ways due to the soft character of LC 
matrix. Different responses are obtained if NPs have substantial impact on LC orientational 
ordering. For example, let us assume that a spherical NP enforces radial distribution of 7 at the 
NP-LC interface. Therefore, it locally acts as a topological defect represented by m=1, see 
Figure Ic. If the total topological charge of the system equals zero (this is the case if far from 
NP liquid crystalline ordering is preferently homogeneously aligned along a single direction), 
then LC matrix must respond by creating antidefects (for example one with m=-lor two with 
m=-1/2) to compensate the effective charge introduced by the nanoparticle. 
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Furthermore, in 2D systems it has been demonstrated that the local surface curvature can 
have strong impact on position and even number of TDs (Kralj, et al., 2011, Vitelli, Turner, 
2004). In particular, the Gaussian curvature G plays a significant role. It can be expressed with 
local principal radia R, and R, of curvature asG =1/(R,R,). It has been shown that G has 
similar impact on topological defects as electrical field on electrical charges. Namely, regions 
with positive (negative) value of G attract TDs with negative (positive) charge m. Furthermore, 
surfaces exhibiting regions with G<0 and G>0 have potential to generate pairs {defect, antide- 
fect} (Vitelli, Nelson, 2006, Vitelli, Turner, 2004, Jesenek, et al., 2013). Therefore, interactions 
between TDs and NPs within a LC on curved surfaces (e.g. in LC films or within biological 
membranes possessing) can give rise to diverse complex structures. 

One way to obtain different properties of soft-matrices is to follow examples 1n music. 
In the latter case on uses few different music notes and their structural organization defines a 
melody. In the case of soft-matrices role of notes can be played by topological charges. They 
can be placed on different positions by playing with surface local curvature and/or appropriate 
nanoparticles. In the following we use a simple model to illustrate key mechanisms via which 
positioning of “notes” in orientationally ordered soft materials can be controlled. 


Minimal Model 


A minimal mesoscopic theory describing nematic LC orientational ordering in 2D curved 
surface embedded in 3D coordinate system have been used. Points on a surface are parametrized 
by parameters (u,v). Details are given in (Kralj, et al., 2011). In reality such LC configurations 
could be encountered in thin LC films, biological membranes exhibiting in-plane orientational 
order or other analogous systems. 

We set that a local surface orientation is determined by its outer unit normal v . A local 
average orientation of a LC rod-like molecule is described byn . The corresponding tensor 
order parameter is defined as 


O=|(7@n-i, @n,). (1) 
where J stands for the Identity tensor, 7 xn, =V, and the scalar order parameter |_ measures 
the degree of ordering. 

The corresponding orientational free energy density fcan be expressed as a Landau-type 
expansion in terms of Q taking into account symmetry allowed terms. In the lowest order 
approximation it is given by (Kralj, et al., 2011) 


f =-ATrQ’ + B(TrOY +kIV, Q| (2) 

















Here A, B, k are positive material constants and ||. marks the surface gradient operator. 

In elastically undistorted bulk system the equilibrium degree of ordering in orientationally 
ordered phase is given byl . = ,/A/(2B) . The competition between the condensation (the first 
two terms in Eq.(2)) and the elastic (the last term in Eq.(2)) free energy contributions defines 
the order parameter correlation length of the system, which is roughly given by x ~Vk/ A. 
Equilibrium LC configurations are obtained via numerical mimization of the free energy for a 
given 2D geometry (Kralj, et al., 2011). 


Results of Theoretical Case Study 
In the following we demonstrate strong impact of curvature and also simple NPs on 
position and number of topological defects of charge m= +1/2 which play in our analogy the 


role of “notes”. We refer to m=1/2 and m=-1/2 TDs as defects and antidefects, respectively. For 
demonstration reason we use as simple as possible model capable to display variety of different 
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structures which have potential to strongly influence resulting physical properties. For this 
purpose we consider thin nematic films covering objects exhibiting spherical topology. 
Furthermore, we restrict to axially symmetric shapes. For better visualization we also consider 
geometries where typical linear dimension R of a shape is larger but comparable to x ~ Vk / A. 
The latter length roughly determines the size of the core of defect, where degree of ordering is 
locally strongly distorted. Note that the core of a defect roughly corresponds to the “body 
shape” of “notes”. Namely, in most situations the core structures (at least degree of ordering) 
are relatively weakly influenced by details of surrounding LC ordering (Kralj, Virga, 2001). 

In Figure 2a we show positioning of TDs on an ideal sphere. Note that due to topological 
reasons such geometry enforces total topological charge equal to m,=2. Exactly at a defect 
origin the nematic ordering is melted (i.e., | =0), making them clearly visible in our simula- 
tions. Relative position of defects and their number and core structure is better visualized by 
plotting degree of ordering in the (u,v) plane in Figure 2b and Figure 2c. In the former we plot 
| (u,v) and in the latter the director field. Here wu €[0,2p[ and ve[0,p[ describe the azi- 

muthal and the zenithal angle of the spherical coordinate system. Figure 2c clearly re- 
veals that four TDs bearing m=1/2 are formed. Note that in this case the Gaussian curvature is 
constant over the sphere. For this reason position of TDs is dominated by mutual repulsion 
among defects. Consequently, TDs occupy vertices of tetrahedral in order to maximize their 
mutual separation (Vitelli, Nelson, 2006). 
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Figure 2: Topological defects in orientational order within a thin LC film cover- 
ing the surface of a sphere. a) Degree of ordering superimposed on the 
sphere. Defects are localized at regions exhibiting black color. The bar at the 


right hand determines the value! /I ..b) Two dimensional variation of | /I . 
in the (u,v) plane, where u and v describe the azimuthal and zenithal angle of 
the spherical coordinate system. c) The corresponding director field distribution 
in the (u,V) plane. 
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In Figures 3 we demonstrate impact of curvature on positions of defects. For this 
reason we consider prolate and oblate ellipsoids. In these cases positive curvature exhibits 
maximum at the poles and equatorial ring, respectively. With increasing prolateness pairs of 
defects are approaching poles (Figure 3a). On the contrary, increasing oblateness the defects 
are approaching the equatorial line (Figure 3b). Note that if we enforce two spherical NPs at 
the poles of a sphere which impose radial director field at their boundaries director field would 
be aligned along meridians. In this case TDs are absent. Namely, such NPs effectively act as 
TDs each bearing topological charge m=1, imposing the total topological charge m, =2 to the 
system. Consequently, in this case “real” TDs are absent. 

Structures shown in Figures 4 exhibit spatially dependent Gaussian curvature and possess 
both regions with G>0 and G<0. Negative value of G 1s established at the necks of structures. In 
Figures a) to d) we progressively increase the average mean curvature of structures. In Figure 
4a and 4b there are four defects as in the spherical case (Figure 2a). However, due to spatially 
nonhomogeneous dependence of G the defects are pushed towards the poles of objects. For 
strong enough local curvatures new pairs defect-antidefect can appear. In Figure 4c we show 
the configuration bearing eight TDs. Three defects are positioned at each bead and at each neck 
one antidefect resides. In Figure 4d we have 12 TDs. 





Figure 3: Prolate and oblate ellipsoids. Defects with positive value of m tend to 
move towards regions exhibiting maximal positive value of the Gaus- 
sian curvature. 
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Figure 4: Shapes exhibiting spherical topology enforcing m,,,=2. Average cur- 
vature of structures gradually increases from a) to d). Configura- 
tions in a) and b) possess 4 defects with m=1/2. Within structures 
in c) and d) curvature is strong enough to generate additional pairs 
(defect, antidefect). 


Summing-up 


Music is composed of a finite set of different music notes. Importantly, by varying 
structure formed by these notes qualitatively different music styles could be formed. Music 
is a typical example where a structure defines a property. This concept could be followed in 
education when nano-sciences are the subyect of lecture. In presented contribution relevant 
basic mechanisms in nematic liquid crystals, which are typical representatives of soft materials 
had been presented. 

If thin 2D nematic LC films have been considered the degree of orientational ordering 
with the tensor order parameter could be described. For demonstrative purpose could be used 
the spherical topology, which by itself enforces the total topological charge m,, =2, which 
showed that by manipulating a local geometry (i.e. the Gaussian curvature G) a number and 
position of topological defects bearing topological charge m =+1/2 could be controlled. Here 
TDs in LCs mimic notes in music. Research demonstrated that TDs can be efficiently targeted 
to different positions within a LC matrix by playing with local curvature. In particular defects 
(TDs characterized by m=+1/2) are attracted to regions exhibiting the maximal positive 
value of G. On the contrary antidefects (TDs characterized by m = —1/2 ) tend to be assembled 
at regions exhibiting minimal negative value of G. Furthermore, by imposing regions exhibiting 
both positive and negative values of G arbitrary number of additional TDs could be introduced 
provided that a LC region is large enough. Positions and number of TDs could be in addition 
manipulated by appropriate NPs. For example, if a nanoparticle imposes a radial-type director 
profile in its close neighborhood it effectively affects surrounding LC medium like a topological 
defect bearing topological charge m=1. First of all, such NP can strongly influence positions of 
neighboring TDs. Next, due to the topological charge conservation these NPs can change 
number of TDs. Therefore, by exploiting geometry and NPs different spatial patterns of “notes” 
(i.e. TDs) can be formed. Different patterns could in general enable different physical responses. 
In LCs, which are famous for their tunability and responsiveness to local stimuli, one could 
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realize rich variety of patterns of TDs exhibiting qualitatively different symmetries. Note that 
several examples in nature demonstrate that symmetry strongly influences physical properties. 
For example, diamond, carbon nanotubes, graphene and graphite are all composed of Carbon 
atoms and they only differ only in its structural spatial (1.e. symmetry) arrangement. 

As a final Remarque, objects of our study, named also LC shells (Vitelli, Nelson, 2006, 
(Kralj, et al., 2011, Nelson, 2002), could be exploited as basic unit elements for future soft col- 
loidal crystals. Here different arrangement of colloids within the crystal would result in differ- 
ent physical properties in a similar way as different packing of atoms results in different crystals 
made of real atoms. Namely, Nelson (Nelson, 2002) proposed to disperse appropriately coated 
LC shells within isotropic media including appropriate nano-sized binders. The latter would 
anchor to defect sites and in this manner organize shells into regular crystal superstructures. 
This research topic is still in its infancy. At this stage only various defect structures in relatively 
simple geometries (spherical and elliptical) have been analyzed. There is a need to find simple 
ways to control sensitively the valence of LC shells and in particular to develop strategies to 
assemble them 1n crystal structures of desired symmetry. This would allow tailoring specific op- 
tical dispersion relations or other physical property of interest and make possibilities to explain 
students very complicated structures in a user friendly way. 
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